A Few Lessons from pQCD Analysis at Low Energies!: 

D.V. Shirkov 
Abstract 

Motivated by the recent 4-loop analysis of the JLab data on Bjorken Sum Rule, where 
the pQCD series seems to blow up at \Q\ < 1.5 GeV, Us > 0.33, we overview the general 
origin of the divergency of common perturbation expansion over powers of a small coupling 
parameter in QFT and consider in detail the blowing-up phenomenon and accuracy of finite 
sums for simple alternating and non-alternating examples of divergent series. 



Introduction 



a; 



It is known since the mid-XX that the main computational tool of quantum theory, 
the perturbation expansion ^fc«'^Cfc(...) over powers of the small coupling parameter 
a , is not a convergent one; expansion coefficients grow factorially Ck ^ k \ . The reason 
is that every quantum amplitude (matrix element) C{a, ...) is not a regular function 
of a at the origin a = . 

Practically, the finite sum J2k such, a series could blow up at iV ~ 1/a . 
To illustrate, take a formal divergent series 



n>l 



Its finite sum 



f[k]i.9) = /ni fn 



nig' 



(2) 



according to the Poincare estimate can 
approximate an expanded function F with 
accuracy Akf{g) = f{g) - f[k]ig) ~ fk ■ 




Fig. 1 Values of fk terms at g=0.25. 



Thus, the finite sum can provide us with the best possible accuracy [AKf{g)]opt 
fx at an optimal number of terms 



k = K r^l/g. 



(3) 



The very existence of this lower limit of possible accuracy is an exact antithesis to the 
case of convergence series : any attempt to increase the number of terms above K leads 
to the lower accuracy. At (7 < 1 this can happen for rather small K values. 

In the above formal example (1), at g=0.25, with K=4, /4(0.25) = 0.5625 
and A4F(0.25) = fi = 3/32 this lower limit of accuracy is about 16.7%. For 
/5(0.25) = 0.6798, it is slightly worse - 17.2 %. 



^The text of contribution to Proceedings of "Intern. Workshop on e+e Collisions from cj) to J/^*" 
(PHIPSIll), Novosibirsk Sept 2011; to be published in Nucl.Phys.(Proc.Suppl.) 
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2 Divergent Series and their Summation 



2.1 Explicit Illustrations 

Consider the integral 

POO 

A{g) = j^ / e-^^-(^/^)^^rfx; g>0. (4) 
Jo 

Expanding integrand in g and changing the order of integration and summation 
one arrives at alternating divergent series 



e""' x^" dx; Ao = l. (5) 

n>0 

The 7T, — )■ oo limit for An coefficients can be estimated by the steepest descent method: 



A"' ~ / e"^(^) dx, /(x) = 4 Inx ; with result A° 

Jo n 



Here, the divergent series was obtained by formal manipulation with the finite 
expression. The finite sums a[n]{g) = g Ai — ■ ■ ■ ± An {—gY of alternating series ([5]) 
can be compared with exact values of the functioij^ A{g) = 1 — a{g) . For results of 
comparison see Figj^l There, we show that starting from g = 0.25 the 0^4] curve passes 
farther from the exact one than the a[3] curve. 




A practical example of alternating divergent series gives the beta-function of the 
gif^ model. In the late 70s its expression 

C(^7) = i/-f^7' + 16.27/ -135.8/ 

^Expressible via the particular Bessel function A{g) — e^l"^^ (tt^)"^/^ Kii^i^jlg) with known an- 
alytic properties. It is analytic in the whole g complex plane (cut along the negative real semi-axis) 
with essential ^ e^^/^ singularity at the origin ; for details, see Sect. 2.2 in paper [2]. 
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known up to the N^LO term was used [2] as a starting point for the whole function 
/^Ms(fi') restoration. In the reconstruction procedure (based also on asymptotic expres- 
sion [4J for 13'^^ ) the Borel representation supplemented by conformal transformation 
was involved. The resultinai closed formula 



9 



dx J 



(6) 



can be used for next coefficient estimation. Later on, the next N'^LO term was calcu- 
lated [3 E] /Ss = 1420.6 via Feynman diagrams. Comparing it with the prediction (6) 
(3^^ = 1409.6 gives the accuracy within 1 % ! 



Another model integral 

Ci9) = ^ 



^ ""^^ dx 



g^Ck ; Ck — Ak 



r(2fc+i/2) 
4'=r(fc+i) 



r(fc) 



(7) 



produces non-alternating asymptotic power series with the same coefficients. As far as 
this integral is also expressible in terms of Bessel functions (see Ref.[2], page 482), one 
has exact expression for the coefficients and can compare the finite sum approximations 
C[ni(fi') = Ci + ■ ■ ■ + CnQ"^ with exact values of C{g) — 1 = c{g) - see Fig.|3l 
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Fi gure 3: The cjj,] approximants for the function c{g) — C{g) — 1 . 

It is clear that the 2-term approximant (lower thin curve) is good up only to g = 
0.15 — 0.20 and the 3-term one (upper thin curve) up to (7 ~ 0.33 while the 4-term 
sum (upper broken curve) starts to deviate from C{g) (red thick curve) at (? ~ 0.27! 

The model (7) is more instructive for our case motivated by the fresh signal from the 
perturbative Quantum Chromodynamics (pQCD) in the low-energy domain. There, 
the 4-loop analysis of rather precise JLab data on polarized Bjorken Sum Rule revealed 
[7] that the non-alternating series for the pQCD correction (eq.(3) in [7J) 



0.3183 as + 0.3631 a? + 0.6520 + 1.804 at 



^with 63(2;) being the cubical polynomial in w{x) , the conformal variable. 
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does blow up at ~ 0.35 . It is noteworthy that the coefficient ratios here (1.1, 1.8, 2.8) 
are close to the factorial ones (1,2,3). 
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Figure 4: The as -dependence of the relative PT contributions to the Bjorken ampUtude at the 
4- loop level - based on Eq.(9). 



Indeed, as it is shown on FigJH the 4-loop term 



a; 



is close to the 3-loop one 



in the interval 0.3 ^ ^ 0.4 , while it approaches the 2-loop term at > 0.4 . We 
marked the first region as a "yellow zone" and the second as a "red" one. Roughly, 
this corresponds to the rule if ~ l/a^ , Eq-(3) with K = i — \ . 

Two other illustrations on Fig. 5 a,b, also taken from paper f7] demonstrate the lack 
of progress in the 4-loop description - marked by black hatching (SW-NE direction)- 
with respect to the 3-loop one (red hatching in the NW-SE direction). 
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Fig. 5a The QCD perturbation analysis of 
the Bjorken form-factor confronted with 
JLab data in three- and four-loop orders. 
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Fig. 5b Instability of the HT coefficient /i4 
fitting of the Jlab data, as in Fig. 5a 



2.2 Asymptotic series and essential singularity 

Turn to the origin of the non- convergent asymptotic series (AS) like in Eqs.(l),(6),(8). 
Usually, it is related with the essential singularity at the origin a = that is a common 
property (in the theories of Big Systems) of the objects representable via Functional 
or Path Integral. This is the case for Turbulence, Classic and Quantum Statistics 
and Quantum Fields. Numerous examples are well known : the er^l^ dependence of 
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the energy gap in BCS and Bogoliubov theories of SuperConductivity : for tunnehng 
probabihty in quantum mechanics. In the theory of Quantum fields (QFT) it was first 
discussed for QED by Dyson just 60 years ago [S] and soon after that implemented by 
Bogoliubov f^; (for the QCD case, the same method was used in the so-called APT 
approach - see below Section 3). 

Mathematically, the essential singularity origin is connected with the small param- 
eter g (or a ) attached to some nonlinear structure. In the quantum case, this is 
interaction term. Generally, a certain AS corresponds to a set of various functions. 
Hence, in physics. 

The Asymptotic Series "summation" is an Art. 

This motto really implies that for the adequate AS summation one should involve 
some additional arguments, like in the Eq. (6) example above. 

2.3 Higher PT terms for e+e — > hadrons 

As far as this meeting is devoted mainly to the electron-positron collider physics, 
turn to the inclusive e^e" — t- hadrons process. Two functions, the cross-section ratio 
R{s) = 1 + r(s) , and the Adler function Di^Q"^) = 1 + d{Q) are in use there. Table 
1 presents the short summary of the PT terms relative contribution in the 'moderate 
energy' interval below m-r . 

Table 1. Relative contributions of 1-, ... 4-loop terms in e+e~ — hadrons 



Function 


Scalc/Gcv 


PT terms (in %) 


the loop number I — > 


1 


2 


3 


4 


r(s) 


1 


65 


19 


55 ?!? 


-39 ?!? 


r(s) 


1.78 


73 


13 


24 ?!? 


-10 ?!? 


d(Q) 


1 


56 


17 


11 ! 


16 ?!? 


d(Q) 


1.78 


75 


14 


6 


5 ! 



In the upper two lines, for r(s) , one can see the literally terrible effect of the tt^ 
terms on the higher £ = 3,4 contributions. This issue was resolved in the 80s[Tni ITT] . 
The net result is that in the annihilation channel, the s-channel, one should use some 
special QCD coupling a{s) instead of as(Q^) . See below, eq.(9) and Fig.6a. 

Concerning the higher contributions, ^^=3,4 , to the Adler function one observes the 
picture analogou^ to the one illustrated by Fig. 4. 

3 Analytic Perturbation Theory 
3.1 A Few Words about APT 

Analytic Perturbation Theory (APT) in QCD, is the closed theoretical scheme devisecl^ 
in the mid-90s [12] without Landau singularities and additional parameters. It stems 

^ with due account of the QCD common coupUng values as{\ GcV) = 0.55 and as{mT-) — 0.35 . 
^See also review papers [T^ \TM [T5] . 
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from the imperatives of RG-invariance, Q^-analyticity, compatibility with hnear inte- 
gral (like, the Fourier) transformations and essentially incorporates non-perturbative 
e"^/"" (algebraic in Q^j^ structures. 

Instead of the power PT set as{Q'^),Cis{Q^)'^, c^siQ'^)^, ■ ■ ■ one has a non-power 
APT expansion set {AkiQ"^) } k = 1,2, . . . with all AkiQ"^) regular in the IR region. 
Accordingly, for the s-channel, there is another IR- regular set ak{s) . The first functions 
^i(Q^) = ctan 5 and ai{s) = a{s) at the one-loop case look rather simple 

«an(Q') = i^(g2/A2) + (A2 _ Q2) ^ = ^TCt^^ ^^^^y (9) 

Both are presented on Fig. 6a together with common , singular at Q = A = 400 MeV . 
Their regular LE behavior corresponds quantitatively to results of lattice simulation 
(see Fig.6b) down to Q ~ 500 MeV. 



Asymptotic fit of 
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Fig.6a Analytic QCD couplings aa„(Q) and pig gb The lattice based on thrcc-gluon 

a(s) in comparison with common cxg ■ vertex 



As it can be seen from Fig. 7 the APT+Higher Twist (HT) description of the JLab 
data looks quite satisfactory down to 350-400 MeV , that is to the Aqcd scale ! 

We omit here technical details of the APT+HT analysis of paper [7]. Some of 
them can be seen in the last right columns of Table 2. There, higher PT and APT 
contributions to couple of sum rules are summarized. 



Table 2. Relative contributions (in %) of l-,2-,3- and 4-loop terms 



Process 


Scale/Gev 


PT(in %) 


APT* 


the loop number = 


1 


2 


3 


4 


1 


2 


3 


Bjorken SR 


t 


1 


35 


20 


19 ! 


26 ?!? 


80 


19 


1 


Bjorken SR 


t 


1.78 


56 


21 


13 


11 ! 


80 


19 


1 


GLS SumRule 


t 


1.78 


65 


24 


11 ! 




75 


21 


4 


Incl. T-decay 


s 


1.78 


51 


27 


14 


7 ! 


88 


11 


1 



* The 4-loop APT contributions are negligible everywhere. 



^For the deep connection between the a -non-perturbativity and the Q^-analyticity, see Ref.[3] 
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Figure 7: Reasonable fit of the JLab data by APT supplemented by Higher Twist (HT) terms (the 
upper green dotted curve) down to 350 McV . Three lower curves describe the standard PT fits. 

Invitation for Work (instead of Conclusion) 

A number of topics is in order: 

• Devising methods of AS summation, (including integral and conformal tricks), 

• Devising Generating Function for HT terms in QCD 

• either generalizing the minimal APT, 

• Toy models for the 4-loop term predicting for other processes Pi 

• Set of analytic couplings a\ , each being adequate to a given process Pi ? 

• Generating HT function for the each Pi ? 
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